Light Vector Mesons at Finite Baryon Density 



T. Hatsuda 

Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 



Abstract 

We summarize the current theoretical and experimental status of the spectral 
changes of vector mesons (p, a;, (p) at finite baryon density. Various approaches 
including QCD sum rules, effective theory of hadrons and bag models show 
decreasing of the vector meson masses in nuclear matter. Possibility to detect 
the mass shift through lepton pairs m. ^ — A, p — A and A — A reactions are 
also discussed. 

I. INTRODUCTION 

At high temperature (T) and density (p) , hadronic matter is expected to undergo a phase 
transition to the quark-gluon plasma. The order parameter characterizing the transition is 
the chiral quark condensate (gg), the absolute value of which decreases as (T,p) increases. 
Numerical simulations of quantum chromodynamics (QCD) on the lattice are actively pur- 
sued to determine the precise nature of the transition at finite T |I| and various model 
calculations have been done to look for the observable signature of the phase transition. 

In this talk, I will concentrate on one of the interesting critical phenomena associated 
with the QCD phase transition, namely the spectral change of hadrons, in particular the 
mass shift of light vector-mesons (p, uj and 0) in nuclear matter at zero T. The vector mesons 
are unique in the sense that they decay into lepton pairs {e^e~ and fi~^^~) which can be 
detected experimentally without much disturbance by complicated hadronic interactions. 

In section 2, I will review the current knowledge of the quark condensate in medium. 
In section 3, various approaches to calculate the vector meson masses in nuclear matter 
are summarized. Experimental possibilities to detect the spectral change are discussed in 
section 4. Concluding remarks are given in section 5. 



II. QUARK CONDENSATES IN NUCLEAR MATTER 

The medium modification of the quark condensate has been calculated since then by 
various methods (lattice QCD, chiral perturbation theory, Nambu-Jona-Lasinio model etc). 
See an overview [0. By these studies, it turned out that there is one noticeable difference 
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between the behavior of {qq) at finite T (with p = 0) and that at finite p (with T = 0): In 
the former case, the significant change of the condensate can be seen only near the critical 
point T ~ Tc 0. On the other hand, in the latter case, 0(30%) change of {qq) could be seen 
even in normal nuclear-matter density. This observation is based on the following formula 
in the fermi-gas approximation (independent particle approximation) |^ 

{uu) 4S^jv fPF d^p Mn , . 

{uu)o'^ flmlJ (27r)3E(p) " ^ > 

Here mM{mT,) is the nucleon (pion) mass, is the pion decay constant, Et^at = (45±10)MeV 
is the TiN sigma term, and E{p) = \Jp^ + M^. (■) and denote the expectation value 
in nuclear matter and that in the vacuum respectively. The integration for the nucleon 
momentum p should be taken from to the fermi momentum p^. At normal nuclear matter 
density {p = Po = 0.17/fm^), the above formula gives (34±8)% reduction of the chiral con- 
densate from the vacuum value. In Fig.l, {uu) / {uu)o as well as the strangeness condensate 
{ss)/{ss)o are shown in the linear density approximation where the uncertainty of E^rAr 
is considered. Estimates taking into account the fermi motion and the nuclear correlatons 
show that these corrections at p = po are less than the above uncertainty @]. 

Unfortunately, the condensate itself is not a direct observable and one has to look for 
physical quantities which are measurable and simultaneously sensitive to the change of the 
condensate. The masses of light vector-mesons are the leading candidates of such quantities. 



III. VECTOR MESONS IN NUCLEAR MATTER 

Let's consider p, u and mesons propagating inside the nuclear matter. Adopting the 
same fermi-gas approximation with (|l|) and taking the vector meson at rest (q = 0), one 
can generally write the mass-squared shift as 

where fvN{p) denotes the vector-meson (V) - nucleon (N) forward scattering amplitude 
in the relativistic normalization, and my(my) denotes the vector meson mass in nuclear 
matter (vacuum). Here, we took spin-isospin average for the nucleon states in fvN- If one 
can calculate fvN{p) reasonably well in the range < p < Pp = 270 MeV (or 1709 MeV < 
^/s < 1726 MeV in terms of the V — N invariant mass), one can predict the mass shift. 
Unfortunately, this is a difficult task: fvN{p) is not a constant in the above range since 
there are at least two s-channel resonances A^(1710), A^(1720) in the above interval and two 
nearby resonances iV(1700) and A (1700). They all couple to the p — N system 0] and give 
variation of fvN{p) as a function of p in principle. From this reason, one should develop 
other methods to estimate 6my without refering to the detailed form of /yiv(p). We will 
briefiy review two of such approaches in the following subsections, namely the QCD sum 
rules and effective theories of hadron. 
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A. QCD sum rules 



The QCD sum rules (QSR) for vector mesons in nuclear matter were first developed by 
Hatsuda and Lee 0. In their approach, one starts with the retarded current correlation 
function in nuclear matter, 

U,,{u,q) = I J d^xe''>^RJ^{x)JM) , (3) 

where g'^ = {uj, q) and RJ^(x) J;y(0) = 9{x^)[J^{x), Ju{0)] with the source currents defined 
as J^''^ = ^{u'jfj^u =F dlfid) (— (+) is for the p°(u;)-meson) and = sj^s. Although there 
are two independent invariants in medium (transverse and longitudinal polarization), they 
coincide in the limit q — and reduce to n^^/(— 3c<j^) = 11. 11 satisfies the following 
dispersion relation, 

1 r°° Imllfii) 
ReUiuj^) = -P / du^ ^ ^ i + (subtraction). (4) 

TV Jo — UJ^ 

In QSR, the spectral density Imll is modeled with several phenomenological parameters, 
while Rell is calculated using the operator product expansion (OPE). The phenomenolog- 
ical parameters are then extracted by matching the left and right hand side of (|^) in the 
asymptotic region —>■ —oo. The density dependence in the OPE side is solely determined 
by the density dependent condensates which are evaluated from low energy theorems or from 
the parton distribution of the nucleon ^j. 

In the medium, we have three kinds of structure in the spectral density: the resonance 
poles, the continuum and the Landau damping contribution. For q — >■ 0, the last contribution 
is calculable exactly and behaves like a pole at cj^ = ll,^. In total, the hadronic spectral 
function looks as 

87rlmn(^z > 0-) = 6{u')p,, + F*6{u^ - ml*) + (1 + ^)9{u' - S*) = PhadXu^), (5) 



with = 27r^p/yp|, + ~ 2n'^p/MN. rriy, F* and Sq are the three phenomenological 
parameters in nuclear matter to be determined by the sum rules. 

Matching the OPE side and the phenomenological side via the dispersion relation in 
the asymptotic region u"^ — oo, we can relate the resonance parameters to the density 
dependent condensates. There are two major procedures for this matching, namely the Borel 
sum rules (BSR) [|l^ and the finite energy sum rules (FESR) [ill], which can be summarized 
as 

roc 

/ dsW{s)[pHadXs)-p^,,{s)]=0, (6) 


W{s) = e{So - s) (FESR), e-'/^' ( BSR). (7) 

Here the spectral function phadX^) stands for eq.(||). Popsi^) a hypothetical imaginary 
part of n obtained from OPE. 

To make quantitative analyses of spectral parameters, the stability analysis based on 
the Borel transform is more suitable than FESR. Since the Borel mass M is a fictitious 
parameter introduced in the sum rule, the physical quantities should be insensitive to the 
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change of M within a Borel interval Mmin < M < Mmax! namely the principle of minimum 
sensitivity (PMS) is used. One can accomplish this insensitivity by choosing suitably 
at given density. In Fig. 2, the Borel curves for the p{uj) meson for three different values 
of baryon density are shown with Sq chosen to make the Borel curve as flat as possible in 
the interval 0.41GeV^ < < l.SOGeVl The upper (lower) bound of the Borel interval 
is determined so that the power (continuum) correction after the Borel tranform does not 
exceed 30 % of the lowest order term in OPE. 

By making a linear fit of the result, one obtains [0] 




1 - (0.16 ±0.06) — , (8) 



Po 

(0.15 ±0.05) — , (9) 
Po 

— = 1 - (0.24 ±0.07) — , (10) 

Po 



and 



^ = 1 - (0.15 ±0.05) y —, (11) 

Po 

where y is the OZI breaking parameter in QCD defined as y = 2{ss)j\f/ {uu + dd) a? with (■) 
being the nucleon matrix element, y takes the value 0.1 — 0.2 [Q]. The decrease in eqs. (^pAf) 
is dictated by the density dependent condensates (gg), {{qq)^) and {q'-y^D^q). The errors 
in the above formulas are originating from the uncertainties of the density dependence of 
the these condensates. The contribution of the quark-gluon mixed operator with twist 4, |^ 
which may possibly weaken the mass shift, is neglected in the above. Shown in Fig. 3 is the 
mass shift given in eqs. (|§,|n|) with possible theoretical uncertainties. 

Some sophistications of the QSR analyses by Hatsuda and Lee have been done later by 
several authors. 

(i) Asakawa and Ko have introduced a more realistic spectral function than (|^) by taking 
into account the width of the p-meson and the effect of the collisional broadening due to 
the n — N — A — p dynamics ||12|. By doing the similar QSR analysis as above, they 
found that the negative mass shift persists even in this realistic case. The width of the rho 
meson in their calculation decreases as density increases, which implies that the phase space 
suppression from the p ^ 27r process overcomes the collisional broadening at finite density. 
Further examination of this interplay between the mass shift and the collisional broadening 
is important in relation to the future experiments. Also, finite temperature generalization 
of the Asakawa-Ko's calculation should be done. 

(ii) Monte Cairo based error analysis was applied to the Borel sum rule by Jin and Leinweber 
||13| instead of the Borel stability or PMS analysis employed in [Q. They found rn*^^/mp^^ = 
1 — (0.22 ± 0.08)(p/po) and m'^/m^ = 1 — (0.01 ±0.01)(p/po), which are consistent with eqs. 
(§,|ll|) within the error bars. 

(iii) Koike analysed an effective scattering amplitude fvN defined as 5rn?^ = fvN ■ P using 
the QSR in the vacuum |jl^. Although his original calculation predicting fvN > is in 
error as was pointed out in ref. H,|r5l) revised calculation gives a consistent result with eqs. 



|,|Tl|) within the error bars [0. Note here that fvN does not have direct relation to the 



scattering length at zero momentum /yAr(O). 
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B. Effective theories 



There have been many attempts so far to calculate the spectral change of the vector 
mesons using effective theories of QCD. The first attempt by Chin |T^ using the quantum 



hadrodynamics (QHD) shows increasing u;-meson mass in medium due to a process analogous 
to the Compton scattering; 

uj + N ~^uj + N. (12) 

For the p-meson, similar but more sophisticated calculations taking into account A-resonance 
and in-medium pion show a slight increase of the p-meson mass In these calculations, 
only the polarization of the Fermi sea (the particle-hole excitations) was considered. Also 
their predictions are different from the general assertion by Brown and Rho claiming that 



all the hadron masses except for pion should decrease [18 . 

On the other hand, Saito, Maruyama and Soutome, and Kurasawa and Suzuki have 
realized that the mass of the a;-meson is affected substantially by the vacuum polarization 
of the nucleon in medium 

a; ^ N*N* u, (13) 

where N* is the nucleon in nuclear matter which has smaller effective mass than that in the 
vacuum. They show that the vacuum polarization dominates over the Fermi-sea polarization 
in QHD and leads decreasing vector meson mass. This conclusion was later confirmed by 



several groups |20] and was generalized for the p and mesons |21]. Jaminon and Ripka has 



also reached a similar conclusion by using a model of vector mesons coupled to constituent 



quarks |^ 



Saito and Thomas have examined a rather different but comprehensive model (bag model 
combined with QHD) and found decreasing vector-meson masses [^]; ^pi^/'^p.o; ~ 1 — 
0.09(p/po)- The spectral shift of the quarks inside the bag induced by the existence of 
nuclear medium plays a key role in this approach. 

Basic idea common in the approaches predicting the decreasing mass may be summarized 
as follows. In nuclear matter, scalar (a) and vector [u) mean-fields are induced by the 
nucleon sources. These mean-fields give back-reactions to the nucleon propagation in nuclear 
matter and modify its self-energy. This is an origin of the effective nucleon mass < 
in the relativistic models for nuclear matter. The same mean-fields should also affect the 
propagation of vector mesons in nuclear medium. In QSR, the quark condensates act on the 
quark propagator as density dependent mean-fields. In QHD, the coupling of the mean-field 
with the vector mesons are taken into account through the short distant nucleon loop with 
the effective mass M^. In the bag-model, the mean fields outside the bag acts on quarks 
confined in the bag and change their energy spectrum. 

Let us show here that one can understand the negative mass shift of the vector mesons in 
a simple and intuitive way in the context of QHD. More quantitative discussion will be given 
in the later section. After renormalizing infinities in the vacuum loop, the density-dependent 
part of the Dirac-sea polarization to the vector-meson propagator is approximately written 
as 

- l^-TT T = 2 2 ^ (14) 

Z ^q^ — rriy — Zniy 
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where Z being the finite wave-function renormahzation constant in medium. The pole 
position is thus obtained as rriy = ^/Zmy- Because of the current conservation, only the 
wave function part of the propagator is modified in medium. Since the effective mass of the 
nucleon decreases in medium {M^/M^ < 1)5 physical vector mesons have more probability 
to be in virtual baryon— anti-baryon pairs compared to that in the vacuum. This means 
Z <1, which leads to m*y/mv = Z <1 |20|j2l|. 

M^/Mn < 1 ^ Z <1 m*y/mv = Z <1 . (15) 



IV. EXPERIMENTS 

How one can detect the spectral change of vector mesons in experiments? One of the 
promising ideas is to use heavy nuclei and produce vector mesons inj — A or p — A reactions. 
Suppose one could create a vector meson at the center of a heavy nucleous. (It does not 
matter whether it is created at the nuclear surface or at the center as far as the produced 
vector mesons run through the nucleous before the hadronic decay). It is easy to see that the 
number of lepton pairs decaying inside the nucleous Nin{l^l~) and that outside the nucleous 
Nout{l^l~) are related as 

Nout{in-)^ e-r-« ' ^ ^ 

where Ttot denotes the total width of vector mesons ((1.3fm)~^, (23fm)~^ and (45fm)~^ for 
p, uj and 0, respectively) and R being the nuclear radius. Eq.(|l^) shows that even the (j) 
meson has considerable fraction of Nm/Nout if the target nucleous is big enough. 

There exist already some proposals to look for the mass shift of vector mesons in nuclear 



medium p4[. One is by Shimizu et al. who propose an experiment to create p and ui in 
heavy nuclei using coherent photon - nucleus reaction and subsequently detect the lepton 
pairs from p and u. Enyo et al. propose to create meson in heavy nuclei using the proton- 
nucleus reaction and to measure kaon pairs as well as the lepton pairs. By doing this, one 
can study not only the mass shift but also the change of the leptonic vs hadronic branching 
ratio r = r(0 — > e"'"e~)/r(0 — > K~^K~). Since is very close to 2mx in the vacuum, any 
modification of the 0-mass or the K-mass changes the ratio r substantially as a function of 
mass number of the target nucleous. Similar kinds of experiments are also planned at GSI. 

There are also on-going heavy ion experiments at SPS (CERN) and AGS (BNL) where 
high density matter is likely to be formed. In particular, CERES/NA45 and HELIOS-3 at 
CERN reported enhancement of the lepton pairs below the p resonance pB|,|2B|, which may 



not be explained by the conventional sources of lepton pairs. If this phenomena is real, low 
mass enhancement of the lepton pair spectrum expected by the mass shift of the vector 
mesons could be a possible explanation . In nuclear collisions at higher energies (RHIC 



and LHC), hot hadronic matter or possibly the quark-gluon-plasma with low baryon density 
are expected to be formed. In such cases, double 0-peak structure proposed by Asakawa 
and Ko as well as the spectral change of p, ui and scalar mesons ||29| will be a distinct 



signal of the chiral restoration in QCD. 
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V. CONCLUDING REMARKS 



The spectral change of the elementary excitations in medium is an exciting new possibility 
in QCD. By studying such phenomenon, one can learn the structure of the hadrons and the 
QCD ground state at finite (T, p) simultaneously. Theoretical approaches such as the QCD 
sum rules and the hadronic effective theories predict that the light vector mesons (p, uj and 
0) arc sensitive to the partial restoration of chiral symmetry in hot/dense medium. These 
mesons are good probes experimentally, since they decay into lepton pairs which penetrate 
the hadronic medium without loosing much information. Thus, the lepton pair spectroscopy 
in QCD will tell us a lot about the detailed structure of the hot/dense matter, which is quite 
similar to the soft-mode spectroscopy by the photon and neutron scattering experiments in 
solid state physics. 
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Figure Captions 

Fig.l The light quark condensates in N=Z nuclear matter in the linear density approxima- 
tion. Theoretical uncertainty of the tiN sigma term is taken into account. We take 
y = 0.12 for the OZI breaking parameter, where y = 2{ss) n/ {uu + dd)^ with 
being the nucleon matrix element. 

Fig.2 Borcl curve for the p{uj) meson mass. Solid, dashed and dash-dotted lines correspond 
to p/ pq = 0, 1.0 and 2.0 respectively. Sq{p) determined by the Borel stability method at 
each density is also shown in GeV^ unit. The Borel window is chosen to be 0.41GeV^ < 
< 1.30GeVl 

Fig.3 Masses of p, uj and mesons in nuclear matter predicted in the QCD sum rules. The 
hatched region shows theoretical uncertainty. 
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